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ABSTRACT 


Dual-Frequency Polarized Scatterometer (DFPSCAT) is a 
new system utilizing Doppler beam sharpening (DBS) 
technology for azimuthal resolution enhancement. 
Considering the DBS technology is inapplicable for the 
middle areas of the swath, a theoretical framework of 
deconvolution signal processing is proposed to improve 
resolution. A deconvolution method of the nonlocally 
centralized sparse representation (NCSR) is adopted to 
verify its feasibility, and simulation results show that the 
deconvolution method have obviously better resolution 
enhancement and higher recovery accuracy than these of the 
classical scatterometer image reconstruction (SIR) method. 


Index Terms—DFPSCAT, resolution enhancement, 
deconvolution method 


1. INTRODUCTION 


A scatterometer is a radar device that collects normalized 
radar backscatter ( o ) measurements over the earth. 
Scatterometer is originally intended to measure near-surface 
wind vector over the ocean. Due to its frequently global 
coverage, it is also a desirable candidate for land and ice 
applications, in which the improvement of spatial resolution 
is the key issue. The pencil-beam scatterometers are real- 
aperture systems and employ the range-gate discrimination 
to enhance the slant resolution; nevertheless, the resolution 


in the azimuth direction is limited by the antenna beamwidth. 


The limitation has been overcome partially with the advent 
of scatterometer image reconstruction (SIR) algorithm, 
which has been successfully utilized to detect polar sea ice 
extent, track icebergs and other applications [1]. However, 
the resolution is still too coarse for many geophysical 
applications, such as snow water equivalent (SWE) 
estimation. In this situation, spatial resolution is desired to 
be 1-5km [2]. 

To achieve the desired spatial resolution, the Doppler 
Beam Sharpening (DBS) technology together with the pulse 


compression technology are adopted for the dual frequency 
polarized scatterometer (DFPSCAT) onboard the Water 
Cycle Observation Mission (WCOM) satellite [3]. 
DFPSCAT is a pencil-beam scatterometer employing a 1.5m 
dish antenna operating at 9.6 and 17GHz for SWE 
measurement and frozen thorn detection. However, the DBS 
is inapplicable for the middle areas of the swath where the 
azimuth resolution is limited by the antenna beamwidth and 
it is about a few tens kilometers. Respecting the low 
resolution of the middle swath of DFPSCAT, the 
deconvolution signal processing method is explored to settle 
the problem. 

The paper establishes a scatterometer measurement 
model, and proposes a new surface o° grid definition 
convenient for solving the model through some 
deconvolution methods. A deconvolution algorithm of 
nonlocally centralized sparse representation (NCSR) is 
utilized to realize the resolution enhancement. Some 
simulations are carried out to verify the feasibility. The 
results show that the proposed framework has obvious 
resolution improvement than that of the well-founded SIR 
processing. 


2. MEASUREMENT MODEL FOR DFPSCAT 


2.1. Measurement model 
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Fig. 1 DFPSCAT measurement geometry 


Figure 1 shows the DFPSCAT measurement geometry. The 
antenna beam is conically scanning around the nadir axis 
with the predesigned rotating speed at a constant incidence 
angle 6. The measurement power P, from a transmitted 
pulse equals the sum of the backscattered signal power P, 
and the noise power P,, which is introduced by thermal 
noise in the receiver, radiometric noise and speckle. The 
power reflected from the surface is 
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where P, is the transmitted power, 4 is the wavelength, L is 
the system loss, G is the antenna gain, o° denotes the true 
normalized radar backscatter of the earth’s surface, R is the 
slant range, and A is the area illuminated by the antenna, v 
represents a defined coordinate system. 

Using the radar equation, a weighted spatial average of 
is derived as 
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h represents the spatial response function (SRF) defined by 
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Generally, P, is estimated through subtracting P,, which 


is measured by the noise-only channel from P, Besides, the 
calibration error resulting from the uncertainty of radar 


system parameters also introduces errors between ō° and 


dv (1) 
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where X = dv , is the calibration factor, and 


measured ó’. The total errors n=G°-G° is characterized by 
N(0, K7e°) , where 


the scatterometer noise model, i.e.n 


K, is the normalized standard deviation of ô’. Then, we 
obtain the following scatterometer measurement model. 


6 =0 +n=| hod +n (4) 
2.2. Surface o° grid 


In the above measurement model, the observed °’ is made 
with a sampled aperture approach in which the 
measurements are spatially filtered surface o° data sampled 
over a two-dimensional (2-D) grid. 

Typically, the 2-D grid is defined as a uniform grid 
takes the along track and cross track direction of satellite as 
coordinate axis [4]. Fig. 2 depicts the typical grid definition, 
where r is a distance from nadir point to center of SRF, and 
the azimuth angle g denotes the instantaneous scan position 
of the antenna beam. Fig. 3 illustrates three SRF A at g =77°, 
83°, 90°, respectively. Note that h varied with each rotation 
location. 6° is not equal to a convolution operation of a 
fixed h and o°, and thus some existing deconvolution 
method are not available for the case. 
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Fig. 2 Typical grid definition 
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Fig. 3 An illustration of SRF for several antenna rotation 
angle. 


To settle the problem, we proposed a new grid 
definition. The new grid is composed of iso-range lines and 
the iso-azimuth lines that the beam boresight line projected 
on the surface, as shown in Fig. 4. Fig. 5 depicted the 
corresponding h with the same locations as Fig.3, and h is 
almost same for varied scan angle. In the case, the measured 
data ° is equivalent to the regularly sampling of surface 
o° convolved with a fixed A, that is 

& =h@o°+n (5) 
where ® denotes the convolution symbol. 
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Fig. 4 A new grid definition 
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Fig. 5 An illustration of SRF in the new grid at azimuth 
angle as Fig. 3 


3. DECONVOLUTION METHOD 


Typically, the DFPSCAT backscatter return from each pulse 
of each beam is partitioned into several “slices” using range 
gate processing. Assuming each slice is statistically 
independent and has its own slice spatial response function. 
Based on the above proposed 2-D grid, the slice spatial 
response function is fixed for measurements. Then, the 
measured G° corresponds to each slice can be written as 

& =h @o° +n, k=1,2++m (6) 
where m is the slice number and h‘ represents the spatial 
response function for kth slice,. Using multiple overlapping 
slice measurements, the surface o° field is reconstructed 
through the deconvolution method. 

Due to the ill-posed nature of deconvolution problem, 
the regularization-based techniques are used to regularize 
solution spaces. Since the classic regularization models like 
Tikhonov regularization and the total variant (TV) 
regularization result in over-smooth solutions, the recent 
sparsity-based regularization used in image restoration is the 
optimal alternative. In the paper, the nonlocally centralized 
sparse representation (NCSR) model proposed by Dong is 
adopted to solve Eq. (6) [5]. Specifically, the high- 


resolution &° is obtained through 


6° =arg min fros eof -> lo, -2 L, (7) 
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where 2, is defined as the weighted average of a,. 

Let o? = Po” be a patch extracted at position i (i =1, 
2, ..., N), where P; is the patch extraction operator and N is 
the number of pixels in the o°. Given a dictionary D, the 
patch o° can be sparsely encoded over D and result in a 
sparse coding vector æ, . Once the coding vectors of all o° 


patches are obtained, the whole o° can be reconstructed by 


o°=Deall (PTR ij >” P" Da, (8) 
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where œ is the concatenation of all @,. More detailed 
explanations on NCSR can be found in [5]. 


4. PERFORMANCE ANALYSIS 


To verify the feasibility of NCSR deconvolution method for 
scatterometer resolution enhancement, some simulations are 
carried out. In the simulation, the slice response function 
contours are approximately 1km in the range direction and 
16km in the azimuth direction. The azimuth-sampling 
interval is set as 1km. In the paper, K, is assumed as 0.5dB. 
At first, a synthetic o” field reconstruction was taken 


as an example for performance investigation. For the o° 
field over the proposed 2-D grid, the left squares are 
1kmx1km in size and spaced 2km apart, the right square is 
2km width, and the width of the narrow lines and broad ones 
correspond to Ikm and 2km, respectively. For comparison, 
the measurement G° and the reconstructed &° after SIR 
algorithm are shown. 

Figure 6 shows the simulation results for DFPSCAT 
resolution enhancement. As expected, the reconstucted &° 
after SIR is excessively smooth and seriously spread out. 
The result of the proposed deconvolution algorithm has the 
most obvious resolution enhancement; nevertheless, it tends 
to more noise amplification compared with the others. 

To future investigate the processing accuracy, the 
horizontal cut of values extracted from Fig. 6 across the top 
row of spots is analyzed, as shown in Fig. 7. It is apparent 
that the proposed deconvolution method can reconstruct the 
truth signal with highest accuracy. 


(b) 


(d) 


Fig. 6 Simulation results for DFPSCAT resolution 
enhancement (a) synthetic o° field (b) measurement (c) 
SIR reconstruction (d) NCSR deconvolution 
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Fig. 7 Plot of o° values extracted from top row of Fig. 6 


Second, a real o° field reconstruction is implemented 


for further validation. The o° field is obtained by a SWE 
forward model of DMRT-QCA [6] with SWE distribution 
data centered at 49°N, 115°W with areas of 2°x2° for 
February 29 in 2004, where the SWE distribution data is 
acquired from Snow Data Assimilation System [7]. The 
reconstruction results are displayed in Fig. 8. It is noted that 
the NCSR deconvolution methods discriminates the edge 
areas clearly and preserves more detail information. Overall, 
the deconvolution method produces the more obvious 
resolution enhancement and higher recovery accuracy 
compared with the classical SIR processing. 


(a) (b) 


dB 
-10 
-15 
-20 
-25 
(d) 


Fig. 8 Real o° field reconstruction (a) real o° field (b) 
measurement (c) SIR reconstruction (d) NCSR 
deconvolution 
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5. CONCLUSION 


This paper introduces a deconvolution theoretical 


framework for scatterometer o° resolution enhancement, in 
which the slice measurement model of scatteometer is 
established based on a new defined 2-D surface o° grid. A 
deconvolution algorithm of nonlocally centralized sparse 
representation (NCSR) is utilized for the DFPSCAT middle 
swath to verify its feasibility. Some simulation experiments 
are conducted and the results show that the reconstructed 


o” after the deconvolution method produces obvious better 
resolution enhancement and higher recovery accuracy 
compared with the classical SIR algorithm. 
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